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The  work  described  in  this  volume  was  performed  under  Contract  No. 
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This  report  describes  research  conducted  on  the  nitrogen  and  phos¬ 
phorus  budget  of  the  Lake  Conway  ecosystem  and  is  part  of  the  Large- 
Scale  Operations  Management  Test  (LSOMT)  of  the  Aquatic  Plant  Control 
Research  Program  (APCRP)  at  the  WES.  A  series  of  WES  reports  document 
the  findings  of  this  research.  Report  1  presents  baseline  information, 
and  Reports  2  and  3  contain  data  collected  the  first  and  second  years, 
respectively,  after  stocking  the  lake  with  the  white  amur.  The  work 
reported  herein  is  the  fourth  volume  (of  seven)  in  Report  2. 

The  work  was  performed  and  this  volume  was  written  by  Dr.  Eldon  C. 
Blancher ,  II,  and  Mr.  Charles  R.  Fellows  of  the  University  of  Florida, 
Gainesville,  Fla. 

The  work  was  monitored  by  the  WES  Environmental  Laboratory  (EL) , 
Dr,  John  Harrison,  Chief.  The  study  was  under  the  general  supervision 
of  Mr.  Bob  0.  Benn,  Chief,  Environmental  Systems  Division  (ESD) ,  EL. 

Mr.  J.  Lewis  Decell  was  Manager,  APCRP.  Principal  investigators  at  WES 
for  the  study  were  Mr.  Eugene  G.  Buglewicz  and  Dr.  Andrew  Miller  of  the 
ESD,  EL. 

Commanders  and  Directors  of  WES  during  the  conduct  of  the  study 
and  preparation  of  the  report  were  COL  John  L.  Cannon,  CE,  COL  Nelson  P. 
Conover,  CE,  and  COL  Tilford  C.  Creel,  CE.  Technical  Director  was 
Mr.  F.  R.  Brown. 

This  report  should  be  cited  as  follows: 

Blancher,  E.  C.,  II,  and  Fellows,  C.  R.  1982. 
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LARGE-SCALE  OPERATIONS  MANAGEMENT  TEST  OF  USE  OF  THE  WHITE 


AMUR  FOR  CONTROL  OF  PROBLEM  AQUATIC  PLANTS 
FIRST  YEAR  POSTSTOCKING  RESULTS 

Nitrogen  and  Phosphorus  Dynamics  of  the 
Lake  Conway  Ecosystem:  Loading  Budgets 
and  a  Dynamic  Hydrologic  Phosphorus  Model 


PART  I :  INTRODUCTION 

1.  The  process  of  eutrophication  or  enrichment  in  natural  waters 
is  caused  primarily  by  the  increased  availability  of  plant  nutrients. 

In  most  lake  systems,  the  elements  usually  responsible  for  the  accelera¬ 
tion  of  this  process  are  nitrogen  and  phosphorus.  The  high  concentra¬ 
tions  of  these  nutrients  in  eutrophic  lakes  may  reflect  either  high  load¬ 
ings  from  external  sources  or  alteration  of  nutrient  cycles  within  a 
lake.  Such  enrichment  modifies  the  character  of  a  lake  by  increasing 
primary  productivity  and  thus  leads  to  problems  associated  with  increased 
growth  of  vascular  plants  or  algae.  In  order  to  fully  understand  these 
processes,  it  is  necessary  to  develop  nutrient  budgets  from  which  the 
supply,  losses,  and  transformations  of  nutrients  can  be  ascertained. 

2.  The  work  presented  herein  is  part  of  the  Large-Scale  Opera¬ 
tions  Management  Test  (LSOMT)  being  conducted  by  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  Vicksburg,  Miss.,  to  test  the  effi¬ 
cacy  of  the  white  amur  ( Ctenophavyngodon  idella)  as  an  aquatic  plant  con¬ 
trol  agent  in  Lake  Conway,  Florida.  In  addition  to  an  assessment  of  the 
impact  of  this  fish  on  the  aquatic  macrophyte  populations,  data  on  physi¬ 
cal  and  chemical  parameters,  as  well  as  on  nontarget  biota,  have  been 
collected  to  identify  additional  effects  of  the  white  amur.  The  analyses 
in  this  report  serve  as  a  base  from  which  the  overall  impact  of  this 
fish  on  the  Lake  Conway  ecosystem  can  be  ascertained. 

3.  The  objectives  of  this  investigation  were  to: 

a.  Assess  external  nitrogen  and  phosphorus  loadings  to  the 
Lake  Conway  ecosystem. 
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b.  Determine  if  and  when  nutrient  limitation  exists  in  the 
lake  and,  if  so,  which  nutrient  is  limiting. 

_c.  Develop  a  dynamic  hydrologic — materials  model  for  the  Lake 
Conway  system  to  determine  if  the  budgets  are  reasonably 
accurate. 

<i.  Assess  the  importance  of  internal  nutrient  loadings. 
Objectives  a  and  t>  were  met  and  presented  in  a  preliminary  report 
(Blancher  et  al.  1978).  The  work  reported  herein  presents  the  revised 
nutrient  budgets  and  the  dynamic  model  that  were  subsequently  developed 
to  complete  the  objectives  c  and  d. 

4.  A  number  of  additional  studies  recently  have  been  published  on 
the  Lake  Conway  LSOMT.  These  studies  include  a  general  description  of 
the  Lake  Conway  area  (Theriot  1977)  and  studies  of  its  fishes  (Guillory, 
Land,  and  Gasaway  1977),  aquatic  macrophytes  (Nall,  Mahler,  and  Schardt 

1977)  (Nall,  Schardt,  and  Burkhalter  1978),  plankton  and  benthos  (Fox 
et  al.  1977;  Conley  et  al.  1978),  and  water  quality  (Sawicki  1977),  as 
well  as  system  modeling  (Ewel  and  Fontaine  1977;  Fontaine  and  Ewel  1978), 
zooplankton  and  trophic  state  (Blancher  1979) ,  seepage  flows  (Fellows 

1978) ,  and  a  study  of  the  nitrogen  cycle  (Sompongese  1978).  The  reader 
is  referred  to  these  articles  for  additional  discussion  of  these  topics. 
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PART  II:  METHODS  AND  MATERIALS 


Hydrology 

5.  Lake  and  watershed  areas  were  calculated  from  U.  S.  Geological 
Survey  (USGS)  topographical  maps  with  an  electronic  planimeter  (Hewlett 
Packard  9810A-9864A) .  Extent  of  catchment  area  was  determined  from 
topographic  features  and  by  consulting  area  tax  assessment  maps  at  the 
Orange  County  Public  Works  Office  for  street  drainage  patterns.  Land 
use  and  development  patterns  for  the  watershed  were  estimated  from  recent 
(1975)  aerial  photographs  supplied  by  county  engineers  (Orange  County 
Public  Works,  Orlando,  Fla.). 

6.  A  bathymetric  map  was  produced  using  both  depth  profiles  ob¬ 
tained  with  a  recording  depth  meter  and  an  existing  bathymetric  map.* 

The  area  of  each  metre-depth  interval  was  determined  with  an  electronic 
planimeter,  and  lake  volumes  were  calculated  by  summing  the  volumes  at 
each  1-m  interval.  Mean  depths  of  the  pools  were  computed  by  dividing 
their  volumes  by  their  respective  surface  areas. 

7.  Monthly  rainfall  and  evaporation  data  were  obtained  from  U.  S. 
Weather  Bureau  climatological  data  reports  (National  Oceanographic  and 
Atmospheric  Agency  (NOAA)  1972-1978).  Rainfall  data  reported  by  the 
Orlando  weather  station  at  McCoy  Air  Base  were  used  because  of  the  close 
proximity  of  the  airport  to  the  lake.  Averages  of  evaporation  pan  data 
from  the  Lisbon  and  Lake  Alfred  stations  (within  a  50-mile  (80.5-km) 
radius  of  Lake  Conway)  were  assumed  to  be  representative  of  the  Lake 
Conway  area.  Coefficients  for  estimating  evaporation  from  the  lake 
surface  were  those  determined  for  Lake  Okeechobee  by  Kohler  (1954) . 

8.  The  volume  of  direct  runoff  was  calculated  by  the  Department 
of  Agriculture  Soil  Conservation  Service  (SCS)  method  and  the  rational 
method  (Chow  1964) .  The  SCS  method  is  based  on  an  empirical  model  that 
relates  the  amount  of  direct  runoff  to  land  use  and  hydrologic  soil 


*  Schardt,  J.,  Unpublished  bathymetric  map  of  Lake  Conway,  Department 
of  Natural  Resources,  Tallahassee,  Fla. 
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characteristics.  A  weighted  curve  number  (CN)  is  developed  from  tab- 
ularized  values  for  the  different  soil  groupings  and  the  expected  direct 
runoff  is  obtained  from  the  model.  The  rational  method  is  based  on  the 
assumption  that  a  certain  percentage  of  the  rainfall  for  a  given  area 
will  contribute  to  the  direct  runoff.  By  measuring  the  amount  of  direct 
runoff  during  a  storm  event,  a  runoff  coefficent  can  be  calculated  by 
dividing  the  runoff  by  the  total  volume  of  rainfall  in  the  basin. 

9.  To  calibrate  these  methods,  a  128-ha  area  of  the  drainage 
basin  was  identified  from  topographic  maps  and  confirmed  with  Orange 
county  engineers  for  street  drainage  patterns.  During  an  intense  thun¬ 
derstorm  in  August  1976,  flow  through  an  open  channel  draining  into  Lake 
Conway  from  this  area  was  measured.  Flow  rate  was  measured  with  a 
Gurley  flow  meter  at  15-min  intervals,  and  water  samples  were  collected 
for  turbidity,  specific  conductance,  and  nutrient  determinations.  Rain¬ 
fall  was  simultaneously  collected  and  measured.  Additional  hydrographs 
were  obtained  during  March,  May,  and  June  1978.  Continuous  measurements 
of  flow  rates  during  these  periods  were  performed  using  an  ISCO  bubble 
flow  meter  and  automatic  recorder  in  conjunction  with  a  45-deg  V  notch 
weir  used  as  the  primary  device.  During  each  storm  event  water  samples 
for  nutrient  analysis  were  obtained  with  a  flow-actuated  ISCO  automatic 
sampler . 

10.  Subsurface  seepage  into  Lake  Conway  was  measured  by  Fellows 

(1979),  using  the  methods  of  Lee  (1972,  1977),  and  Lock  and  John  (1978). 

2 

Drums  enclosing  0.255  m  of  lake  bottom  were  pushed  (open  end  down)  10 
to  15  cm  into  the  sediments.  The  closed  ends  were  vented  by  0.9-cm-diam 
plastic  tubing  into  either  a  0.180-  or  3.5-i  bag.  Flow  was  calculated 
by  dividing  the  volume  of  water  collected  by  the  measured  time.  Water 
samples  were  collected  from  the  individual  drum  ends  for  subsequent 
nutrient  analysis.  Drum  collectors  were  placed  on  transects  perpendicu¬ 
lar  to  the  shoreline  since  seepage  flow  decreases  with  increasing  dis¬ 
tance  from  shore  (Lichtler,  Anderson,  and  Joyner  1968).  Flow  as  a 
function  of  distance  from  shore  was  then  plotted  and  integrated  and  ex- 
pressed  as  flow  per  metre  of  shoreline  per  unit  time  (m  /m-day). 


Nutrients 


11.  Dry  fallout  and  rainfall  for  nutrient  analyses  were  collected 
initially  (August  1976-September  1976)  using  large  (4 -l)  nalgene  beakers. 

After  October  1977,  collection  of  these  samples  was  accomplished  using 
an  automatic  wet-dry  precipitation  collector  (Aerochem  Metrics  Inc., 

Miami,  Fla.).  This  collector  consisted  of  two  large  buckets  with  a  servo- 
operated  lid  and  a  precipitation  sensor. 

12.  Beginning  in  June  1976,  monthly  water  samples  were  collected 
for  nutrient  analysis  at  all  stations  in  the  Conway  system  at  a  depth  of 
1  m.  Additional  samples  were  obtained  at  4  and  7  m  at  the  deeper  sta¬ 
tions  in  the  center  of  each  pool.  After  March  1977,  sampling  was 
continued  at  the  deep  stations  only  with  samples  collected  at  each  metre 
interval.  Supplemental  nutrient  data  covering  the  period  January  1976- 
March  1978  were  obtained  from  the  Orange  County  Pollution  Control  Depart¬ 
ment.*  All  water  samples  were  analyzed  for  nitrogen  ^NH*,  NO2,  NO^,  and 
total  Kjeldahl  nitrogen  (TKN)^  and  phosphorus  (ortho  and  total  phosphate) 
using  methods  described  in  Standard  Methods  (American  Public  Health  Asso¬ 
ciation  (APHA)  1975)  and  the  U.  S.  Environmental  Protection  Agency 
Methods  Handbook  (1974) . 

13.  Bioassay  methods  including  algal  assay  procedure  (U.  S.  Envi¬ 
ronmental  Protection  Agency  1971),  alkaline  phosphatase  assay  (Fitzgerald 
and  Nelson  1966)  and  nitrogen  fixation  (Stewart,  Fitzgerald,  and  Burris 
1968)  were  used  to  determine  nutrient  limitation  in  the  Conway  system. 

14.  A  hydrologic  and  phosphorus  model  was  developed  using  a  linear 
systems  model  similar  to  one  described  by  Rich  (1973).  The  model  was 
simulated,  using  the  Continuous  System  Modeling  Program  (CSMP)  (Speck- 
hart  and  Green  1976)  on  the  Northeast  Regional  Data  Center's  Amdahl  460 
computer.  Functions  for  rainfall  and  evaporation  were  derived  using 
data  from  the  National  Weather  Service  (NOAA  1972-1978),  and  lake  height 
function  was  developed  from  the  USGS  (1977-1978).  Data  used  for 

-  0 

*  Personal  Communication,  1978,  Raymond  T.  Kaleel,  Orange  County  Pollu¬ 
tion  Control  Department,  Orlando,  Fla. 
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developing  functions  of  phosphorus  release  from  both  vascular  aquatic 
plants  and  sediments  were  obtained  from  Fontaine  (1978). 
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PART  III:  RESULTS  AND  DISCUSSION 


Development  of  the  Hydrologic  Budget 

15.  Nutrient  inputs  to  a  lacustrine  system  are  closely  related  to 
the  hydrologic  budget  of  the  system.  Therefore,  the  construction  of  an 
accurate  water  budget  is  necessary  for  the  development  of  a  realistic 
materials  budget. 

16.  The  hydrologic  budget  of  a  lake  can  be  expressed  by  the 
following  equation: 


AV=Q.  +P-E-Q  +  Q 
i  o  -  xg 


(1) 


where 

AV  =  change  in  water  storage* 

=  surface  inflows 

P  =  precipitation 

E  =  evaporation 

Q  =  surface  outflows 
vo 

Q  =  groundwater  input  or  output 
§ 

17.  Monthly  precipitation,  evaporation,  and  pan  coefficients  for 

the  1976  water  year  are  presented  in  Table  1,  The  volume  of  direct 

6  3 

precipitation  onto  the  Lake  Conway  surface  was  8.87  x  10  m  .  An  annual 

6  3 

evaporation  volume  of  9.75  x  10  m  was  determined  by  multiplying  the 
monthly  pan  coefficients  (Kohler  1954)  by  the  lake  surface  area 
(7.39  km)  and  the  monthly  pan  evaporation  and  then  summing  the  evapora¬ 
tion  volumes. 

18.  Surface  and  subsurface  inflows  to  the  lakes  were  divided  into 
two  major  components  each:  stream  inflow  and  surface  stormwater  runoff, 
and  subsurface  seepage  and  groundwater  inflows.  Stream  inflow  to  the 
Conway  system  is  primarily  by  small  ditches  that  drain  the  surrounding 
watershed.  Consequently,  inflow  from  these  ditches  was  considered  as  a 


*  All  units  were  10  cubic  metres  per  month. 
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part  of  stormwater  runoff.  Total  runoff  was  calculated  from  stormwater 
runoff  and  seepage.  Stormwater  runoff  is  that  water  which  enters  stream 
channels  associated  with  storms  and  contributes  directly  to  the  runoff 
hydrograph.  Seepage  is  that  portion  of  total  runoff  which  enters  the 
lake  following  infiltration  of  soils  in  the  watershed.  Groundwater  in¬ 
flows  include  those  flows  into  the  lake  from  the  confined  aquifer. 

19.  To  adequately  assess  these  inflows,  it  was  necessary  to  cal¬ 
culate  them  using  appropriate  methods,  combine  them  to  determine  total 
inflows,  and  then  evaluate  how  they  compared  to  expected  total  inflows 
calculated  from  changes  in  lake  stage.  Thus,  a  check  of  the  accuracy  of 
the  measured  inflows  could  be  made.  The  volume  of  stormwater  runoff  can 
vary  considerably  within  a  drainage  basin.  Soil  type,  vegetative  cover, 
and  extent  of  impervious  surfaces  as  well  as  frequency,  intensity,  and 
duration  of  rainfall  play  major  roles  in  determining  the  amount  of  run¬ 
off.  It  is  therefore  necessary  to  obtain  hydrographs  of  stormwater 
flows  for  the  basin  also,  using  standard  empirical  methods  for  deter¬ 
mining  direct  runoff. 

20.  Stormwater  hydrographs  for  the  Conway  basin  were  calculated 
from  data  collected  in  August  1976  and  March,  May,  and  June  1978,  and  are 
summarized  in  Table  2.  An  example  is  given  in  Figure  1,  and  the  remain¬ 
ing  hydrographs  of  these  storms  are  presented  in  Appendix  A.  Measure¬ 
ments  of  stormwater  flows  during  1977  were  not  possible  because  of  the 
paucity  of  rainfall.  Of  the  precipitation  that  fell  during  those  storms, 
3.5  percent  on  an  average  entered  the  lake  as  stormwater  runoff.  Initi¬ 
ation  of  runoff  and  time  interval  to  peak  flow  and  return  to  base  flow 
were  generally  rapid. 

21.  Estimates  of  stormwater  runoff  by  the  SCS  (1975)  method  re¬ 
sulted  in  a  CN  of  56  for  the  Conway  basin.  Using  the  SCS  empirical 
model,  this  means  that  approximately  4.6  cm  of  rainfall  is  necessary 
before  any  surface  runoff  can  be  detected.  For  a  6.35-cm  rainfall 
(2.5  in.),  approximately  0.15  cm  of  runoff,  or  2.4  percent  of  rainfall, 
would  result. 

22.  Since  both  the  rational  and  SCS  methods  gave  comparable  run¬ 
off  coefficients,  the  annual  volume  of  stormwater  runoff  was  calculated 


Figure  1.  Hydrograph  of  stormwater  runoff  on 
21  August  1976  at  the  corner  of  Gatlin  Ave. 
and  Bumby  Dr.,  Orlando,  Fla. 


by  multiplying  the  annual  amount  of  precipitation  in  the  basin  by  0.035, 
the  average  runoff  coefficient  determined  from  the  hydrographs. 

23.  Seepage  flows  into  Lake  Conway  were  measured  for  a  period  of 
12  months  and  plotted  against  the  dates  they  were  obtained  (Figure  2) . 
This  annual  seepage  graph  is  intended  to  reflect  the  total  seepage 
input  to  both  lakes.  The  total  seepage  reflects  not  only  the  steady 
seepage  flows  but  also  some  seepage  that  pulsed  into  these  lakes 
following  rain  events. 
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Figure  2.  Seepage  flows  plotted  against  time  for  the  three  Lake  Conway 
sites  during  1977-1978.  Inset  histogram  shows  monthly  rainfall  values 

during  study  period 


24.  Seepage  measurement  scheduling  was  usually  arranged  2  or  3 
weeks  in  advance,  and  sampling  was  conducted  regardless  of  the  prevail¬ 
ing  weather  conditions.  The  seepage  obtained  during  this  study  is 
believed  to  be  unbiased  with  respect  to  weather  conditions,  and  all 
measurements  were  represented  on  the  annual  seepage  graphs  either  di¬ 
rectly  or  as  an  averaged  value  for  replicate  measurements  of  the  same 
or  consecutive  dates. 

25.  An  annual  seepage  input  value  for  the  East  Pool  site  was 
calculated  by  first  averaging  the  initial  and  final  flow  values  and  mul¬ 
tiplying  this  by  the  9  days  not  included  in  the  356-day  measurement 

3 

period.  This  product  (0.9  m  /m)  was  then  added  to  the  integrated  356-day 
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3  3 

input  (63.9  m  /m  shoreline).  The  annual  East  Pool  input  was  64.8  m  /m 

shoreline. 

26.  The  annual  seepage  values  for  the  remaining  two  Lake  Conway 

sites  were  calculated  by  use  of  a  seepage  quotient.  This  quotient  was 

3 

the  annual  East  Pool  seepage,  64.8  m  /m  shoreline  (integrated  for  1 

year),  divided  by  that  amount  of  East  Pool  seepage  which  occurred  during 

the  time  period  that  each  site  was  measured.  For  the  South  Pool,  the 

integrated  seepage  value  (Figure  2)  and  the  calculated  annual  inputs 
3 

were  53.1  and  98.3  m  /m  shoreline,  respectively.  At  the  West  Pool 

3 

these  values  were  32.3  and  64.2  m  /m  shoreline. 

27.  To  determine  the  average  seepage  value  for  the  five-lake 

3 

system,  the  numerical  average  of  these  three  values,  75.8  m  /m  shoreline 

(integrated  for  1  year),  was  used.  This  average  was  chosen  since  this 

area  was  classified  by  Leighty  et  al.  (1957)  to  be  a  "good"  recharge  area 

for  the  Floridan  Aquifer  and  because  poor  correlations  of  flow  with  land 

use  and  watershed  area  were  made.  From  this  value  and  the  length  of 

shoreline  for  the  five  pools,  the  total  seepage  input  for  the  Lake  Con- 

6  3 

way  system  was  calculated  to  be  2.15  x  10  m  /yr,  which  is  equivalent  to 
30.3  cm  of  lake  stage  per  year. 

28.  The  occurrence  of  a  high  intensity,  short  duration  rain  event 
during  a  routine  sampling  trip  to  Lake  Conway  on  1  December  1977 
afforded  an  opportunity  to  construct  a  seepage  hydrograph.  An  undeter¬ 
mined  amount  of  rain  fell  during  a  10-min  period  immediately  after  the 
routine  sampling  of  the  East  Pool  site.  A  histogram  of  integrated 

flow  versus  time  shows  a  rapid  change  in  flow  and  the  subsequent  decay 

in  rate  to  near  prerain  flow  (Figure  3).  Groundwater  table  heights, 

at  the  times  indicated  (Figure  3) ,  show  a  change  that  corresponds  to 

the  change  in  seepage  flow.  The  integrated  seepage  input  due  only  to 

the  rain  (i.e.,  with  the  initial  seepage  rate  subtracted)  was  13.8  i/m 

_2 

shoreline,  or  7.3  x  10  mm  lake  stage.  To  approximate  the  amount  of 
rainfall,  a  soil  porosity  of  25  percent  was  assumed.  Since  it  was  known 
that  the  maximum  change  in  groundwater  height  was  5.2  cm,  rainfall  was 
then  estimated  to  be  about  1.3  cm.  Seepage  input  from  the  rain  was 
equal  to  0.56  percent  of  what  fell  directly  onto  the  lake  surface. 


TIME  OF  DAY  (PM) 

Figure  3.  Seepage  histogram  and  corresponding  water  table  height  before 
and  after  a  rain  event  at  the  East  Fool  site 
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29.  Estimation  of  the  accuracy  of  the  values  for  the  total  runoff 
components  was  based  on  several  methods  of  determining  total  runoff. 

For  the  first  of  these,  it  was  assumed  that  correlation  of  total  rainfall 
in  the  basin  with  subsequent  rise  in  lake  level  would  provide  an  addi¬ 
tional  estimate  of  total  runoff.  Months  that  contained  well-defined 
storm  events  following  a  dry  period  were  determined  from  the  rainfall 
record,  and  subsequent  changes  in  lake  storage  and  estimates  of  total 
runoff  were  calculated  (Table  3).  From  these  calculations,  it  was  deter¬ 
mined  that  an  average  of  12.7  percent  of  the  precipitation  that  fell  in 
the  watershed  (excluding  the  lake)  entered  the  lake.  Secondly,  informa¬ 
tion  from  the  USGS  (1972-1978)  indicates  that  the  average  annual  outflow 
via  Boggy  Creek  is  17.8  cm.  Considering  that  average  total  rainfall  for 
the  area  is  132.1  cm,  runoff  is  approximately  13.5  percent  of  total  rain¬ 
fall.  Variations  in  rainfall-runoff  for  Boggy  Creek  (Figure  4)  indicate 
that,  when  total  rainfall  was  120  cm,  approximately  12  to  13  percent  of 
the  rainfall  that  fell  in  the  basin  contributed  to  total  runoff. 

30.  By  combining  the  estimated  surface  runoff  and  seepage  input 
and  dividing  by  the  precipitation  in  the  watershed  ((2.24  +  1.62)/ 

3 

46.28  m  ),  a  total  runoff  of  8.3  percent  is  obtained.  Compared  to  the 
average  calculated  total  runoff  above  of  12.7  percent  (Table  3),  this 
value  is  considerably  lower  than  what  would  be  expected  for  this  area. 

31.  An  estimate  of  groundwater  flow  (from  the  confined  aquifer) 
was  determined  by  considering  the  fall  in  lake  stage  during  an  extended 
dry  period  and  correlating  this  value  with  losses  estimated  from  evapora¬ 
tion  alone.  During  three  such  periods,  the  drop  in  lake  levels  could  be 
predicted  by  evaporation  alone,  showing  no  net  inflow  or  outflow  from 
groundwater  and  seepage. 

32.  This  method  does  not  distinguish  groundwater  (confined 
aquifer)  flow  from  subsurface  seepage.  Seepage  of  this  latter  type  was 
not  monitored  during  these  extended  dry  periods  but,  during  this  study, 
was  always  found  to  flow  into  the  system  in  large  amounts  (usually 
>100  l/m  shoreline-day  i.e.,  >0.4  mm/day).  Consequently,  groundwater 
flows  mighc  be  expected  to  be  outflows  from  the  system  since  evaporation 
alone  accounted  for  fall  in  lake  stage.  These  observations,  however,  do 
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Figure  4.  Precipitation-total  runoff  relationship  for 
Boggy  Creek.  Closed  circle  (•)  is  the  datum  point  for 
1974,  a  year  when  the  rainfall  record  was  incomplete, 
possibly  explaining  the  apparent  discrepancy 

not  preclude  the  possibility  of  subsurface  seepage  from  the  system,  es¬ 
pecially  in  the  area  near  the  only  surface  outflow.  These  observations 
do  support  the  hypothesis  that  groundwater  flows  (from  the  confined 
aquifer)  into  the  Conway  system  are  insignificant. 

33.  Surface  outflow  through  the  outlet  during  the  study  occurred 
only  during  August,  September,  and  October  1976.  From  three  measurements 
obtained  during  these  months,  it  was  determined  that  annual  surface 
outflow  was  approximately  1.87  x  10  mJ.  Additional  outflows  from  the 
lakes  might  occur  via  groundwater  recharge,  seepage,  evapotranspiration 
losses,  and  domestic  pumpage.  No  attempt  to  measure  these  flows  was  made, 
and  it  was  assumed  that  these  flows  constituted  the  difference  between 
calculated  inputs  and  outputs. 

34.  From  lake  height  information  for  the  years  1974-1975  and 

1975-1976,  it  was  determined  that  the  net  change  in  storage  in  the  system 

6  3 

for  the  water  year  1976  was  -0.53  x  10  m.  This  amount  was  added  to  the 
outputs  as  additional  seepage,  i.e.  groundwater  recharges,  out  of  the 
lake. 
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35.  Utilizing  all  this  information,  a  hydrologic  budget  for  the 
Conway  system  was  constructed  (Table  4).  It  is  evident  that  the  hydro- 
logic  cycle  of  this  system  is  dominated  by  precipitation-evaporation 
events.  Other  significant  inflows  are  stormwater  runoff  and  seepage, 
with  the  latter  being  the  most  important.  Hydrologic  budgets  calculated 
for  the  individual  pools  in  the  Conway  system  are  presented  in  Table  5. 

Development  of  Dynamic  Hydrologic  Model 

36.  The  calculated  flows  were  used  in  a  dynamic  model  to  predict 
change  in  storage  within  the  lake  system.  An  information  flow  diagram 
for  the  model  was  prepared  (Figure  5) .  The  model  was  simulated  for  the 
period  from  December  1972  through  March  1978  on  a  monthly  basis.  Simu¬ 
lation  of  the  model  from  the  measured  flows  resulted  in  a  large,  unreal¬ 
istic  increase  in  lake  storage  (i.e.,  the  model  was  unstable).  Two  possi¬ 
bilities  could  explain  this  discrepancy:  (a)  either  the  outflows  were 
too  small  or  (b)  the  inflows  were  too  large.  Since  no  reasonable  argu¬ 
ments  could  be  presented  to  justify  additional  outflow  from  the  lake, 

it  was  assumed  that  the  inflows  to  the  lake  were  too  large.  The  control¬ 
ling  factors  in  the  inflow  equation  left  only  two  possibilities  where 
errors  could  exist:  (a)  either  the  seepage  and  runoff  coefficients  were 
too  large  or  (b)  the  basin  area  was  too  large.  Since  earlier  arguments 
indicated  that  the  total  runoff  measured  was  already  too  small,  the 
most  likely  place  of  error  was  the  basin  area.  Even  though  the  area  is 
well  defined  by  drainage  patterns  and  sewers,  the  effective  drainage 
area  is  probably  smaller  because  of  high  infiltration  and  depression 

storage  in  the  area.  Therefore,  the  basin  area  was  determined  to  be 

2 

29.2  percent  less,  and  was  reduced  from  35.6  to  25.2  km  .  As  a  result 
of  this  change,  the  model  provided  more  realistic  predictions. 

37.  Additional  refinement  of  the  model  was  accomplished  by  incor¬ 
porating  a  factor  that  changed  the  runoff  and  seepage  coefficients  to 

reflect  (simulate)  moisture  conditions  in  the  water  table  aquifer.  I 

During  dry  periods  (low  rainfall),  this  factor  decreased  these  coeffi¬ 
cients  and  conversely  increased  them  during  wet  seasons  (high  rainfall). 
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Figure  5.  Flow  diagram  for  the  dynamic  hydrologic — materials 
systems  model.  Qq  =  inflow  volume  at  time  zero,  Qi  =  outflow 
volume  at  time  1,  cq  =  inflow  concentration,  cl  =  average 
lake  concentration,  V  =  lake  volume,  Gy  *  groundwater  vol¬ 
ume,  Qg  =>  groundwater  flow 

Also,  minor  changes  in  the  coefficients  were  necessary  to  obtain  the  best 
simulation.  These  adjustments  resulted  in  a  model  that  simulated  the 
observed  lake  height  accurately  for  the  period  December  1972  to  June 
1977.  After  this  period,  the  model  showed  an  increase  in  lake  height 
that  was  not  observed.  Repeated  ajustments  of  the  important  coefficients 
did  not  improve  the  response  of  the  model.  Since  a  realistic  simulation 
was  desired,  an  additional  drain  (approximately  3  percent  of  inflows) 
was  added  to  reduce  the  levels  to  those  observed  in  the  lake.  The  expla¬ 
nation  that  this  discrepancy  was  due  to  unpredictable  changes  in  outflow 
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characteristics  or  recharge  characteristics  of  the  basin  during  an  ex¬ 
ceedingly  dry  year  is  plausible. 

38.  The  results  of  the  model  are  compared  with  the  actual  changes 
in  lake  height  in  Figure  6.  A  listing  of  the  computer  program  is  in¬ 
cluded  as  Appendix  B.  Differences  between  the  simulated  and  observed 
changes  in  lake  height  are  within  5  percent  on  a  volume  basis. 


Figure  6.  Comparison  of  the  actual  variations  in  lake  height  of  the  Lake 
Conway  system  and  the  simulated  height  from  the  hydrologic  model  for  the 
period  December  1972  (month  0)  through  March  1978  (month  63) 
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PART  IV:  NUTRIENT  LOADING  AND  PHOSPHORUS  DYNAMICS  OF 
THE  LAKE  CONWAY  SYSTEM 


39.  The  trophic  dynamics  of  lake  systems  are  the  integration  of 
incoming  allochthonous  materials  with  the  lake's  biota.  The  first  step 
in  elucidating  factors  affecting  trophic  structure  is  to  quantify  some 
of  these  inputs,  and  this  is  most  effectively  accomplished  by  the  devel¬ 
opment  of  a  nutrient  budget. 

40.  Depending  on  the  type  of  information  desired,  two  approaches 
for  developing  budgets  are  available:  the  net  budget  method  (Johnson 
and  Owen  1971)  and  the  mass  balance  approach  (Burns  1976) .  Net  budgets 
are  based  on  the  concentration  of  the  input  sources  entering  a  system 
and  can  have  an  input  with  a  negative  effect  (i.e.  if  it  has  a  dilution 
effect  on  the  system) .  Net  budgets  are  preferable  when  the  nature  of 
the  input  sources  is  being  investigated.  In  the  case  of  the  mass  budget, 
the  quantity  of  materials  is  considered  independent  of  the  concentration 
and  is  always  positive.  Mass  budgets  are  appropriate  when  the  behavior 
of  the  individual  nutrients  within  the  lake  is  the  prime  consideration 
(Vollenweider  1975).  For  purposes  of  this  investigation,  the  mass  bud¬ 
get  is  the  best  approach. 


Atmospheric  Sources 

41.  Atmospheric  inputs  of  phosphorus  and  nitrogen  originate  from  a 
variety  of  sources  including  fertilizer  mining,  manufacturing,  soils,  and 
combustion  of  various  types  (fuels,  agricultural  burning,  and  forest 
fires)  (Barkdoll,  Overton,  and  Betson  1977).  Eventually  these  materials 
are  removed  from  the  atmosphere  as  dry  fallout  or  scavenged  by  precipita¬ 
tion.  Much  of  the  material  that  falls  on  land  surfaces  can  potentially 
be  resuspended  into  the  atmosphere.  Fallout  onto  water  surfaces  is  es¬ 
sentially  irreversible.  In  these  cases,  open  bodies  of  water  receive  a 
disproportionate  share  of  atmospheric  particulate  matter  (Murphy  1974). 

42.  Based  on  the  rainfall  and  dry  fallout  samples  collected 

2  2 

(Table  6),  annual  loadings  of  0.048  g  P/m  -yr  and  0.36  g  N/m  -yr  for 


l 
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rainfall  and  0.080  g  P/m  -yr  and  0.57  g  N/m  -yr  for  dry  fallout  have 

2 

been  calculated.  The  resulting  bulk  loadings  (0.128  g  P/m  -yr  and 

2  2 
0.93  g  N/m  -yr)  compare  well  with  similar  loadings  of  0.105  g  P/m  -yr 

2 

and  1.0  g  N/m  -yr  for  Gainesville,  Fla.  (Hendry  and  Edgerton  1978, 
unpublished  data,  Department  of  Environmental  Engineering,  University 
of  Florida,  Gainesville,  Fla.). 

Nutrients  From  Stormwater  Runoff  and  Seepage 

43.  The  accumulation  of  particulate  materials  in  urban  areas  and 
its  subsequent  runoff  into  water  bodies  have  been  investigated  recently 
(Fruh  et  al.  i966;  Mattraw  and  Sherwood  1977;  Field,  Curtis,  and  Bowden 
1976).  An  early  study  by  Weibel  (1969)  showed  that  runoff  can  contain 
suspended  solids  at  concentrations  that  exceed  those  in  raw  domestic 
wastewaters.  Seven  probable  sources  of  contaminants  in  urban  runoff  were 
listed  by  Sartor  and  Boyd  (1972):  pavement,  vehicles,  atmosphere,  vege¬ 
tation,  litter,  domestic  and  industrial  spills,  and  antiskid  compounds. 
Transport  of  these  contaminants  into  nearby  waters  often  results  in 
significant  nutrient  loadings  and  subsequent  water  quality  degradation. 

44.  Analysis  of  nitrogen  and  phosphorus  concentrations  in  storm¬ 
water  runoff  from  the  small  catchment  in  the  Lake  Conway  watershed, 
described  in  Part  III,  yielded  average  concentrations  of  3.25  and 

0.39  rng/i  for  nitrogen  and  phosphorus,  respectively  (Table  7).  These 
compare  well  to  published  values  that  range  from  1.93  to  4.45  mg/ l  for 
N  and  from  0.19  to  0.98  mg/i  for  P  (Weibel  1969;  Cowen  and  Lee  1976;  Rast 
and  Lee  1978) .  Multiplying  the  hydrologic  loadings  from  stormwater  run¬ 
off  in  the  watershed  (Table  4)  by  the  measured  concentrations  and  divid¬ 
ing  by  the  surface  area  of  the  lake  resulted  in  areal  loadings  of 

2  2 
0.069  g  P/m  -yr  and  0.84  g  N/m  -yr. 

45.  During  the  period  July  1977  through  March  1978,  seepage 
nutrient  fluxes  were  measured  from  three  to  five  times  at  each  of  the 
sites  in  Lake  Conway.  Table  8  summarizes  the  results  and  shows  the  daily 
fluxes  of  each  nutrient  and  the  seepage  flow  rates  by  site  and  date. 

46.  The  West  and  South  Pool  fluxes  for  August  1977  were  estimated 
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from  incomplete  data  sets.  After  nutrient  flux  patterns  had  been  estab¬ 
lished  for  these  sites,  the  nutrient  fluxes  measured  on  this  date  were 
recalculated  using  graphical  extrapolation  to  include  the  additional  area 
from  the  fourth  seepage  meter  to  the  outermost  meter.  The  estimated 
nutrient  flux  values  probably  still  underestimated  the  true  fluxes. 

47.  Significant  underestimations  of  seepage  nutrient  fluxes  prob¬ 
ably  also  occurred  when  nutrient  fluxes  in  the  meter  furthest  from  shore 
were  relatively  large.  Occurrences  of  this  type  are  indicated  in  Ta¬ 
ble  8.  The  fluxes  of  ammonium  at  all  sites  were  initially  similar,  and 
they  tended  to  decline  after  about  a  2-month  period  of  constant  or  in¬ 
creasing  flux.  The  fluxes  of  ammonium  plus  total  organic  nitrogen  (i.e. 
Kjeldahl  nitrogen  flux)  refined  elevated  for  3  to  3-1/2  months  after 
meter  installation  because  the  flux  organic  nitrogen  was  initially  low 
and  it  increased  as  the  ammonium  flux  declined. 

48.  Total  Kjeldahl  nitrogen  (TKN),  which  accounted  for  most  of  the 
nitrogen  in  the  seepage  water,  remained  relatively  constant  at  each  of 
the  three  sites  in  Lake  Conway  until  after  November  1977.  The  decline 

in  TKN  flux  was  probably  caused  by  lack  of  sedimenting  input  to  the  area 
under  the  seepage  meter,  i.e.,  it  was  an  artifact  of  the  seepage  meter. 

An  average  daily  flux  of  total  nitrogen  was  calculated  for  each  site  in 
the  East  and  the  South  Pools  from  the  data  up  to  and  including  the  Novem¬ 
ber  sampling  date.  Similar  calculations  were  conducted  on  the  ammonium 
and  total  organic  nitrogen  data  from  the  West  Pool  site;  however,  the 
nitrate  plus  nitrite  values  were  not  included  with  this  analysis.  The 
elevated  oxidized  nitrogen  (NC^  and  NO^)  values  were  probably  caused  by 
an  atypically  heavy  fertilizer  application  to  an  adjacent  citrus  grove. 
The  average  of  these  three  site  averages  was  0.62  g  N/m  shoreline-day. 
This  average  daily  flux  was  biased  because  it  was  generated  from  data 
obtained  during  the  months  of  greatest  seepage  flow  and  highest  water 
temperature.  The  effects  of  temperature  on  the  annual  nutrient  flux 
could  not  be  compensated  for  with  the  limited  data  collected  during  this 
study.  The  effects  of  the  seepage  cycle  were  better  known  and  were  used 
to  generate  an  annual  average  daily  flux  value. 

49.  A  seepage  quotient,  similar  to  the  quotients  used  in  the  water 
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budgeting  section,  was  employed  to  calculate  the  annual  average  daily 
flux  of  nitrogen  into  Lake  Conway.  This  ratio  was  generated  from  the 
flow  data  of  the  East  Pool  site,  and  it  consisted  of  the  average  daily 
hydraulic  loading  for  the  year  divided  by  the  average  daily  hydraulic 
loading  for  the  110-day  period  during  which  the  nutrient  flux  data  were 
obtained.  The  product  of  the  average  flux  (0.62  g  N/m  shoreline-day) 
and  the  seepage  quotient  (0.693)  is  an  estimate  of  the  annual  average 
daily  flux;  this  value  was  0.430  g  N/m  shoreline-day.  The  total  annual 
loading  of  nitrogen  from  seepage  was  calculated  by  multiplying  the  latter 

4 

figure  by  the  shoreline  length  of  the  five  pools  (2.95  x  10  m)  and  by 

365  (days  per  year).  The  result,  4640  kg  N/yr,  was  divided  by  the  sur- 

6  2 

face  area  of  the  five  lakes  (7.39  *  10  m  )  to  yield  an  annual  areal 

2 

nitrogen  loading  of  0.628  g  N/m  -yr  for  the  Lake  Conway  system. 

50.  The  phosphorus  loading  to  Lake  Conway  was  calculated  from  the 

mean  of  the  three  site  averages  for  flux  measurements  of  total  phosphorus 

up  to  and  including  the  3  November  1977  date.  The  seepage  quotient 

procedure  was  then  employed  to  minimize  the  effects  of  the  annual  seepage 

flow  cycle.  The  total  phosphorus  seepage  input  to  Lake  Conway  was  180  kg 

2 

P/yr,  or  24.4  mg  P/m  -yr. 

51.  Since  undeveloped  land  in  the  watershed  surrounding  Lake  Con¬ 
way  has  a  high  infiltration  rate,  virtually  no  surface  runoff  should 
originate  from  these  areas.  These  types  of  areas  are  interspersed 
throughout  the  "urban"  area  in  the  drainage  basin.  Consequently,  they 
are  included  in  the  previous  calculations  for  residential  sources. 

52.  Surface  runoff  from  agricultural  areas  is  also  considered 
negligible  for  several  reasons.  Urban  encroachment  into  these  areas  is 
very  evident,  and  their  classification  as  agricultural  is  questionable. 
Few  active  citrus  farms  remain.  Infiltration  in  these  areas  is  also 
great  because  of  the  paucity  of  impermeable  surfaces. 

Additional  Sources  of  Nutrients  to  Lake  Conway 

53.  Other  inputs  of  nutrients  to  lake  systems  occur  in  the  form 
of  internal  loadings:  nutrient  regeneration  from  the  sediments  by 
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chemical  and  biological  processes  (burrowing  insects  and  mixing  by 
fishes);  nutrient  "pumping"  by  vascular  aquatic  plants;  and  nitrogen  fix¬ 
ation  by  blue-green  algae  (and  possibly  by  sediment  bacteria) .  Inputs 
from  the  sediments  and  aquatic  plants  are  not  considered  for  purposes  of 
the  annual  budget,  but  they  will  be  discussed  in  the  phosphorus  modeling 
section.  Detectable  rates  (up  to  30  ng  N/g  dry  weight-hr)  of  nitrogen 
fixation  occurred  during  August  1977  in  Lake  Gatlin  and  in  the  East  Pool 
of  Little  Lake  Conway  associated  with  benthic  algal  mats.  Since  nitrogen 
fixation  only  occurred  at  that  particular  time  and  since  heterocystic 
blue-green  algae  were  not  common  in  the  plankton,  it  was  concluded  that 
fixation  of  nitrogen  represented  an  insignificant  portion  of  the  total 
nitrogen  budget. 


Outputs 

54.  Discharge  of  nutrients  through  the  outlet  of  Lake  Conway 
accounted  for  only  a  small  portion  of  total  nutrient  loss.  It  was  as- 
summed  that  outflow  nutrient  concentrations  were  equal  to  average  surface 
concentrations  (Figure  7),  which  were  0.87  mg/H  for  nitrogen  and 

0.024  mg/Jl  for  phosphorus.  Multiplying  by  the  total  outflow  (Table  4) 

2 

and  converting  to  annual  areal  rates  of  loss  yields  values  of  0.22  g/m  -yr 

2 

for  nitrogen  and  0.006  g/m  -yr  for  phosphorus. 

55.  Losses  by  fish  and  weed  removal,  volatilization,  and  evapora¬ 
tion  were  not  assessed.  However,  losses  via  these  routes  would  probably 
represent  only  a  minor  portion  of  the  total.  Groundwater  recharge  as  a 
nutrient  sink  for  the  Lake  Conway  system  is  considered  to  be  insignifi¬ 
cant,  which  follows  from  the  conclusion  of  insignificant  flow  reached  in 
Part  III.  Even  if  this  conclusion  were  in  error,  nutrients  being  car¬ 
ried  out  via  this  pathway  would  end  up  in  the  sediments  of  the  lake  as 
either  particulates  or  adsorbed  to  clay  and  organic  matter. 

56.  Subtraction  of  nutrient  outputs  from  total  inputs  yields  a 
value  that  is  considered  to  be  a  representative  nutrient  loss  by  sedimen¬ 
tation.  Although  this  loss  does  not  represent  complete  removal  from  the 
system,  sedimentation  does  eliminate  a  significant  portion  of  nutrients 
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from  biological  or  chemical  cycling  within  the  lakes. 

57.  Annual  nitrogen  and  phosphorus  budgets  for  the  Lake  Conway 
system  were  prepared  and  are  presented  in  Table  9.  The  atmosphere  was 
the  major  external  source  of  both  nitrogen  and  phosphorus  to  Lake  Conway, 
accounting  for  37  and  56  percent  of  the  total,  respectively.  Urban  run¬ 
off  was  the  second  most  important  source  of  both  nitrogen  and  phosphorus. 
Monthly  nitrogen  and  phosphorus  loadings  (Table  10)  reflected  differences 
in  the  hydrologic  loadings  to  the  system. 

Phosphorus  Loading  Model 

58.  If  volumetric  flow  rates  into  and  out  of  the  system  are  unequal 
and  time  varying,  and  a  completely  mixed  situation  exists,  the  budget  for 
a  conservative  substance  (not  transformed  chemically  or  biologically  or 
removed  by  sedimentation)  entering  a  lake  would  follow  a  fundamental 

mass  balance  equation  as  follows: 

d(Vc  ) 

-dt”  =  Vi  *  Vo  (2> 


where 

V  =  volume  of  lake  at  time  t* 

c^  =  average  concentration  in  inflow 

=  inflows  to  lake 

Q  =  outflows  to  lake 
o 

cq  =  average  concentration  in  lake 

For  a  nonconservative  substance  such  as  phosphorus,  the  basic  equation 
could  be  rewritten  as: 


d(Vc  ) 

-dt“  =  Vi  -  Vo  -  kcoV  (3) 


* 


All  units  were  10 


6 


cubic  metres  per  month. 
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where  k  is  a  first-order  rate  constant  or  a  term  for  the  loss  of  sub¬ 
stances  other  than  through  the  outlet  (i.e.,  reaction  and/or 
sedimentation) . 

59.  Equation  3  can  be  solved  with  the  restrictive  assumptions 
that  outflow  concentration  is  equal  to  the  mean  lake  concentration 
and  that  sedimentation  is  a  function  of  the  mass  of  material  under  con¬ 
sideration  in  the  lake.  Vollenweider  (1968)  obtained  the  solution: 


w  q  1  +  a  (7/ q  ) 
^s  m  ^s 


(4) 


where 

(m^)  =  mean  lake  concentration  at  steady-state 

=  yearly  areal  loading  of  substance  m 

qg  =  discharge  height,  m 

o  =  sedimentation  coefficient 
m 

~z  =  mean  depth,  m 

Generally,  mean  lake  concentration  is  a  function  of  residence  time,  areal 

loading,  and  sedimentation  characteristics.  Using  this  relationship  for 

phosphorus  in  a  number  of  Swiss  and  North  American  lakes,  Vollenweider 

(1975)  developed  empirical  limits  of  critical  loadings  for  lakes  of 

various  flushing  characteristics.  He  found  that  these  limits  were 
2  2 

0.1  to  0.3  g  P/m  -yr  and  1.0  to  2.0  g  N/m  -yr  in  temperate  lakes. 

2 

Brezonik  and  Shannon  (1971)  found  higher  limits  of  0.28  to  0.49  g  P/m  -yr 

2 

and  2.0  to  3.4  g  N/m  -yr  for  Florida  lakes  using  a  different  approach, 
and  suggested  that  subtropical  lakes  may  be  able  to  assimilate  somewhat 
larger  nutrient  loadings.  Phosphorus  loadings  to  Lake  Conway  fall  below 
Vollenweider ' s  ranges  of  critical  loadings,  whereas  loadings  for  nitrogen 
exceed  the  acceptable  loadings  of  Vollenweider  and  are  in  the  critical 
range  according  to  Brezonik  and  Shannon. 

60.  Two  terms  in  the  Vollenweider  input-output  model  must  be 
carefully  considered  before  utilizing  Vollenweider 's  model:  the  water 
residence  time  and  the  phosphorus  sedimentation  coefficient  . 

Water  residence  time  has  been  recently  shown  to  be  one  of  the  most  im¬ 
portant  factors  in  empirical  loading  models  (Yeasted  and  Morel  1978) . 
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Hence,  critical  to  any  application  of  such  models  is  the  definition  of 

residence  time.  Although  x  appears  to  be  a  straightforward  concept, 

w 

defined  as  the  lake  volume  divided  by  the  outflow,  care  must  be  taken 
on  what  "outflow"  means  in  terms  of  the  models.  As  this  term  represents 
the  flushing  characteristics  of  a  lake,  outflow  (especially  in  cases  of 
lakes  whose  hydrology  is  dominated  by  precipitation  evaporation  events) 
should  be  defined  as  total  outflow  less  the  losses  caused  by  evaporation. 
In  this  way  an  "effective"  residence  time  that  represents  the  true  flush¬ 
ing  characteristics  of  the  system  is  obtained.  All  subsequent  calcula¬ 
tions  were  performed  using  this  "effective"  residence  time. 

61.  The  second  factor  that  must  be  critically  examined  is  the 
phosphorus  sedimentation  coefficient  (equivalent  to  k  in  Equa¬ 

tion  2) .  Realizing  its  importance  in  his  relationship  and  the  reality 
that  it  must  be  calculated  from  the  model  Itself,  Vollenweider  (1975, 
1976)  provided  several  ways  to  calculate  this  coefficient.  The  first 
of  these,  derived  from  observations  on  lakes  in  his  earlier  studies,  is 
as  follows: 


In  cr  *  In  5.5  -  0.85  In  z 
P 


Equation  5  results  in  a  =1.34  for  Lake  Conway.  Use  of  this  value 

P  _ 

to  predict  the  mean  lake  concentration  [P]  from  Equation  3  results  in 
a  value  of  0.034  mg/ 5.  which  is  considerably  higher  than  the  observed 
value  of  0.024  mg/f.  Working  backwards  from  the  observed  concentration, 
a  value  of  1.8  for  a  was  found,  indicating  that  Vollenweider 's  equa- 


a  value  of  1.8  for  a 

P 

tion  can  underestimate 


for  Florida  lakes. 


62.  Another  relationship  which  can  be  used  to  calculate  o 


the  equation  (Vollenweider  1976) : 


KM 


p  *•[?]. 


where 


[P]  =  average  lake  concentration  over  depth  z  ,  mg  P/m 

*  2 
F  (z)  =  the  flux  of  P  through  z  ,  mg  P/m  -yr 


* 


( 


The  major  drawback  to  this  equation  is  that  F^(z)  is  unavailable  for 
Lake  Conway  and  most  lakes  that  Vollenweider  studied. 

63.  Vollenweider  (1976)  then  defined  sedimentation  velocity  as 

and,  by  working  backwards  from  Equation  6  using  a  sedimenta 
tion  coefficient  from  this  model,  he  obtained  an  apparent  settling  veloc 
ity  of  approximately  10  m/yr.  As  Vollenweider  points  out,  the  value  ob¬ 
tained  from  this  model  is  considerably  below  experimentally  obtained 
values  because  it  is  a  net  sedimentation  velocity.  Burns  and  Pashley 
(1974)  measured  actual  settling  velocities  in  Lake  Ontario  and  found 
that  they  ranged  from  -0.4  to  2.0  m/day.  Taking  an  average  of  their  re¬ 
sults  calculated  for  a  year,  a  value  of  240  m/yr  is  obtained. 


A  Dynamic  Model  of  Phosphorus  in  Lake  Conway 


64.  Equation  3  in  the  preceding  section  was  used  in  the  systems 
model  to  simulate  the  change  in  total  phosphorus  concentration  over  time 
(see  Figure  7) .  Although  the  model  was  simulated  for  a  period  of  over 
5  years,  total  phosphorus  concentrations  were  available  only  for  the 
period  January  1976  through  December  1977,  and  for  comparative  purposes 
only  those  results  are  presented  (Figure  7).  The  5-year  simulation, 
however,  was  used  to  test  the  model's  stability. 
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1978  1977 

Figure  7.  A  comparison  of  the  final  simulation  of  the  materials  model 
to  total  phosphorus  concentrations  in  the  lake.  Bars  represent  the 

standard  error  of  the  me  an 

65.  Using  only  external  inputs  and  a  sedimentation  coefficient 
derived  from  Vollenweider 's  (1975)  model,  the  model  yielded  phosphorus 
levels  that  are  representative  of  conditions  in  the  lake  on  an  annual 
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basis,  but  the  simulations  did  not  fit  the  observed  seasonal  trends  of 
total  phosphorus  concentrations.  This  was  expected  because  of  the 
simplifying  assumptions  used  in  Vollenweider ' s  derivation,  i.e.  that  a 
constant  op  can  be  used  to  represent  the  net  sedimentation  coefficient 
for  the  system.  While  the  assumption  of  a  constant  loss  coefficient  may 
be  adequate  for  an  annual  cycle,  it  is  an  oversimplification  for  a 
dynamic  model.  Other  factors  must  be  considered  to  explain  the  seasonal 
events  observed. 

66.  Functions*  were  obtained  for  phosphorus  loadings  from  vascular 
plants  and  the  sediments  from  Fontaine  (1978).  He  estimated  these  func¬ 
tions  from  productivity  measurements  and  sediment  leaching  experiments 
performed  during  1976.  Since  these  functions  represent  maximal  expected 
inputs  from  these  sources,  their  actual  contribution  to  the  lake  is  con¬ 
sidered  as  a  fraction  of  the  tctal.  As  an  initial  estimate  these  func¬ 
tions  were  run  at  half  of  their  potential  loadings.  Reduction  of  the 
loadings  was  continued  until  the  simulated  phosphorus  concentrations  fit 
the  observed  seasonal  trends.  This  occurred  when  the  plants  were  reduced 
to  20  percent,  and  the  sediments  were  reduced  to  1  percent  of  their  re¬ 
spective  maximal  loadings.  In  order  to  obtain  an  optimal  fit,  the  sedi¬ 
mentation  coefficient  from  Vollenweider ' s  model  had  to  be  increased  from 
0.15  to  0.70.  While  some  justification  of  using  a  higher  sedimentation 
coefficient  was  presented  previously  in  the  discussion  of  loading  models, 
it  is  necessary  to  examine  the  relationship  more  closely.  Using  Vollen- 
weider's  relationship  to  compute  a  settling  velocity  for  Lake  Conway 
using  the  sedimentation  coefficient  from  the  dynamic  model,  a  value  of 
44.1  m/yr  is  obtained.  Although  this  is  still  well  below  the  values  re¬ 
ported  by  Burns  and  Pashley  (1974),  it  is  probably  realistic  for  a  lake 
system  with  an  averaged  mean  depth  of  only  5  m. 

67.  As  an  additional  reference  parameter  to  check  the  calculations 
of  the  loading  model,  Vollenweider  (1976)  also  introduced  the  concept  of 
relative  phosphorus  residence  time  ir^  defined  as: 

*  See  Appendix  B,  ***Functions***  of  the  model.  For  time  ■  0.0,  phos¬ 
phorus  (g/m^)  -  0.00102;  for  time  ■  2.0,  phosphorus  ■  0.00108,  etc.; 
time  *  0.1  of  a  year. 
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The  third  term  of  Equation  7  is  developed  from  the  steady-state  model 

and  the  fourth  term  is  a  statistical  approximation.  For  Lake  Conway, 

tt  calculated  from  the  third  term  is  0.053,  but  if  calculated  from  the 
r 

fourth  term  it  is  0.231.  According  to  Vollenweider ,  this  may  indicate 
that  accumulation  of  phosphorus  in  the  sediments  of  Lake  Conway  must 
occur  at  a  greater  rate  than  in  the  lakes  he  used  in  his  study.  This 
corroborates  both  the  results  of  the  dynamic  model  and  the  previous  cal¬ 
culations,  as  well  as  the  findings  of  Brezonik  and  Shannon  (1971).  Shal¬ 
low  subtropical  lakes,  such  as  those  found  in  Florida,  may  be  able  to 
assimilate  more  nutrients  than  would  be  expected  for  temperate  lakes, 
justifying  the  increase  in  the  sedimentation  coefficient. 

68.  The  final  simulation  is  presented  in  Figure  7  along  with  ob¬ 
served  trends  in  total  phosphorus  concentrations.  A  correlation  coef¬ 
ficient  (r  *  0.25)  was  calculated  for  the  two  curves,  indicating  a  weak 
fit  of  the  model  to  the  observed  data.  When  the  first  5  months  of  1976 
are  not  included  in  the  comparison,  a  reasonable  fit  is  observed 

(r  =  0.63).  The  observed  high  phosphorus  concentrations  in  early  1976 
may  have  been  caused  by  a  large  amount  of  phosphorus  being  recycled  in 
the  lake  system  because  of  extensive  herbicide  treatment  of  the  aquatic 
plants  the  preceding  fall. 

69.  Contributions  of  total  phosphorus  of  the  lake  from  the  three 
major  input  sources  are  shown  in  Figure  8  and  in  the  following 
tabulation: 

Year  1976  1977 

External 

All  sources  0.24  0.25 

(Continued) 


*  "w"  refers  to  water.  For  example, 

phosphorus  in  the  lake. 


t  /t  is  the  residence  time  for 
p  w 


% 
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The  relatively  small  contribution  of  phosphorus  from  the  sediments  and 
the  dominance  of  macrophyte  pumping  in  the  model  give  some  insight  into 
the  role  of  these  components,  at  least  for  the  Conway  system.  The  simu¬ 
lation  suggests  that  phosphorus  released  from  the  sediments  during  anoxic 
conditions  in  the  summer  does  not  become  entrained  during  overturn  in 
the  fall.  This  hypothesis  has  been  suggested  previously  by  Fitzgerald 
(1970),  who  observed  that  adsorbtion  of  phosphorus  by  lake  muds  is  rapid 
enough  to  remove  most  of  the  phosphorus  before  it  can  be  used  to  support 
algal  growth.  An  alternate  hypothesis  that  the  amount  of  phosphorus 
from  the  plants  is  too  high  and  the  sediments  are  actually  supporting 
the  observed  concentrations  of  phosphorus  cannot  be  ruled  out  entirely. 
However,  the  mass  balance  of  phosphorus  determined  by  use  of  the  dynamic 
model  does  not  support  the  hypothesis  that  sediment  release  is  an  impor¬ 
tant  source  of  phosphorus  for  Lake  Conway. 


33 


PART  V:  SUMMARY  AND  CONCLUSIONS 


70.  This  report  summarized  and  integrated  some  of  the  physical, 
chemical,  and  biological  data  that  were  collected  from  the  Lake  Conway 
system  during  the  period  April  1976  through  March  1978.  The  purpose  of 
this  work  was  to  clarify  some  of  the  relationships  between  external 
nutrient  loadings  and  internal  nutrient  cycling  in  this  subtropical 
ecosystem. 

71.  Measurements  were  made  of  stormwater  runoff  and  subsurface 
seepage  flows  to  the  lake.  Using  the  collected  data  and  existing  precip¬ 
itation  and  evaporation  data,  a  hydrologic  budget  was  constructed.  The 
budget  was  refined  and  verififed  by  use  of  a  dynamic  systems  model. 

The  hydrologic  model  accurately  predicted  the  change  in  lake  level  for 
the  period  December  1972  through  March  1978. 

72.  Based  on  the  flows  determined  from  the  hydrologic  investiga¬ 
tion,  inputs  of  nutrients  to  Lake  Conway  were  calculated.  The  major 
sources  of  external  nitrogen  and  phosphorus  loadings  were  from  atmos¬ 
pheric  loadings,  urban  runoff,  and  subsurface  seepage.  Estimates  of 

2 

nutrient  inputs  showed  that  both  nitrogen  (2.53  g  N/m  -yr)  and  phosphorus 
2 

(0.224  g  P/m  -yr)  loadings  were  in  the  range  of  "critical"  loading  for 
the  system. 

73.  A  simulation  model  of  the  phosphorus  dynamics  of  Lake  Conway 
was  developed  to  determine  the  relative  importance  of  external  inputs 
and  internal  nutrient  cycling.  From  the  results  of  the  model,  the 
following  conclusions  were  drawn: 

a.  The  sedimentation  of  phosphorus  in  Lake  Conway  occurs  at 
a  rate  higher  than  would  be  predicted  from  Vollenweider ’ s 
analyses  on  temperate  lakes. 

_b.  Nutrient  release  by  submerged  macrophytes  is  an  important 
process  in  the  phosphorus  dynamics  of  Lake  Conway. 

c.  Release  of  phosphorus  by  sediments  is  not  a  significant 
internal  source  of  this  nutrient  for  Lake  Conway. 
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Table  1 

Precipitation,  Evaporation,  and  Pan  Coefficients  for 
Water  Year  1976  (October  1975  through  September  1976) 


Month 

Precipitation* 

cm 

Evaporation** 

cm 

Pan 

Coeff icientt 

October  1975 

12.04 

12.70 

0.76 

November 

1.68 

8.95 

0.71 

December 

1.29 

7.07 

0.83 

January  1976 

0.94 

7.54 

0.77 

February 

2.11 

10.23 

0.69 

March 

4.37 

14.93 

0.73 

April 

5.49 

17.09 

0.84 

May 

26.31 

17.01 

0.82 

June 

25.22 

16.34 

0.85 

July 

17.90 

17.61 

0.91 

August 

8.26 

17.25 

0.91 

September 

14.91 

14.70 

0.85 

Total 

120.52 

161.42 

*  For  Orlando  Weather  Station,  McCoy  Air  Base  station. 

**  Average  for  Lisbon  and  Lake  Alfred  stations, 
t  Kohler  (1954). 


Table  2 

Summary  of  Stormwater  Runoff  Hydrographs  Collected  at  Gatlin  Avenue 
and  Bumby  Drive,  Orlando,  In  1976  and  1978 


Date 

Total 

Flow 

3 

m 

Rainfall 

cm 

Basin 

Area 

ha 

Volume  of 

Rainfall 

in  Basin 

3 

m 

21  August  1976 

986.30 

1.98 

128 

25,344 

23  March  1978 

45.43 

0.18 

43 

774 

4  May  1978 

665.88 

3.68 

43 

15,824 

16  June  1978 

45.75 

0.33 

43 

1,419 

21  June  1978 

4.58 

0.58 

43 

2,494 

Average 

coefficient 


Rain 

Duration 

hr 

Antecedent 
Conditions 
(Time  Since 
Last  Rain) 

Base 

Flow 

m^/ sec 

Runoff 

Coefficient 

0.42 

24  hr 

0.0 

0.039 

1.25 

24  hr 

0.0 

0.059 

2.00 

24  hr 

0.0 

0.042 

0.90 

48  hr 

0.0 

0.032 

0.33 

5  days 

0.0 

0.002 

0.035 


Table  3 

Estimation  of  Total  Precipitation,  Evaporation,  and  Runoff  In 

3  6 

the  Conway  System  (m  x  10  ) 


Date 

Total 
Rainfall 
in  Basin 

A  Volume 
In  Lake 

Precipitation 
on  Lake 

April  1973 

0.686 

+0.112 

0.131 

March  1974 

3.35 

+0.896 

0.638 

May  1975 

Average 

7.63 

+1.299 

1.167 

Evaporation 
From  Lake 

Net 

Precipitation 

Total 

Runoff 

Total  Runoff/ 
Total  Rainfall 
in  Basin 
Expressed  as  a 
Percentage 

0.120 

0.011 

0.100 

15 

0.113 

0.524 

0.371 

11 

0.782 

0.385 

0.914 

12 

12.7 

Table  4 

Annual  Hydrologic  Budget  for  the  Lake  Conway  System 
for  the  Water  Year  1976* 


Volume  (a^  *  10^) 

Sources 

Precipitation  8.87 

Influent  streams  NS 

Stormwater  1.66 

Seepage  in  2.15 

Groundwater  NS** 

(from  confined  aquifer)  _ 

Total  in  12.68 

Sinks 

Evaporation  9.27 

Surface  outflow  1.87 

Seepage  out  NS 

Groundwater  and  pumpage  2.07+ 


Total  out  13.21 

Change  in  storage  -0.53 


*  In  millions  of  cubic  metres. 

**  NS  »  not  significant. 

+  Obtained  by  difference. 


Table  5 

Hydrological  Budget  for  Each  Pool  In  the  Conway  System* 


Pool 

Precipitation 

Stormwater 

Seepage 

In 

Total 

In 

South 

1.64 

0.076 

0.395 

2.11 

Middle 

3.59 

0.136 

0.669 

4.40 

East 

1.52 

0.626 

0.488 

2.63 

West 

1.78 

0.629 

0.433 

2.84 

Gatlin 

0.34 

0.197 

0.164 

0.70 

Groundwater  Total 


Evaporation 

Outflows* ** 

Out 

Out 

1.72 

0.18 

0.21 

2.11 

3.75 

0.20 

0.44 

4.39 

1.59 

0.78 

0.26 

2.63 

1.86 

0.70 

0.28 

2.84 

0.35 

0.28 

0.07 

0.70 

*  In  millions  of  cubic  metres 

**  Obtained  by  difference. 


Table  6 

Loadings  of  Nitrogen  and  Phosphorus  to  Lake  Conway  From  Atmospheric  Sources 


_ Date _ 

18  August  1976* 

19  November-1  December  1977 
1  December-15  December  1977 
15  January-11  February  1978 
19  March-3  May  1978 

3  May-23  May  1978 
23  May-12  June  1978 
12  June-12  July  1978 
12  July-30  July  1978 
Mean  value 


Total  Nitrogen 


Dry 

2 

ig/m  -day 

Wet 

mg/mm  Rainfall 

— 

0.30 

— 

0.25 

— 

0.37 

1.11 

0.17 

1.91 

0.24 

1.58 

0.35 

2.19 

0. 34 

1.32 

0.32 

1.29 

0.30 

1.56 

0.30 

Total  Phosphorus 


Dry 

mg/m^-day 

Wet 

mg/mm  Rainfall 

— 

0.080 

— 

0.030 

— 

0.004 

0.10 

0.007 

0.40 

0.023 

0.18 

0.093 

0.22 

0.030 

0.28 

0.018 

0.15 

0.041 

0.22 

0.036 

*  One  rain  event  only. 


Table  7 

Nitrogen  and  Phosphorus  Concentrations  in  Stormwater 
Runoff  Entering  Lake  Conway 


Number  of 
Samples 

Total 

Nitrogen* 

Total 

Phosphorus* 

Total  Flow 

3 

m 

21  August  1976 

9 

2.3 

0.33 

986.3 

23  March  1978 

27 

4.7 

0.45 

45.5 

4  May  1978 

28 

4.6 

0.47 

665.8 

16  June  1978 

12 

4.0 

0.42 

45.8 

21  June  1978 

8 

5.1 

0.63 

4.6 

Mean  values 

— 

3.25 

0.39 

_ 

*  Concentrations  are  flow-weighted  averages  expressed  as  milligrams  per 
litre. 


Table  8 

Nutrient  Flux  Values  for  the  Three  Study  Sites  In 
Lake  Conway  During  the  1977-78  Study  Period 


Site  and 
Installation 
Date 

Date 

Flow* 

NH*-N** 

4 

East  Pool 

7/16 

196 

0.470 

Some  meters 

11/3 

120 

0.054 

repositioned 

1/6 

103 

0.030 

West  Pool 

8/7 

189 

0. 562++ 

7/27 

9/1 

217 

0.490+ 

11/3 

198 

0.112 

1/6 

209 

0.047 

3/24 

200 

0.015 

South  Pool 

8/7 

283 

0. 525++ 

7/15 

9/1 

368 

0.652 

11/3 

256 

0.359 

1/6 

310 

0.154 

Total 

Organic 

Nitro¬ 

gen** 

N0j-N+ 

N0~-N 

Total 

Phos- 

phorus+ 

Ortho-P+ 

— 

2.57 

11.8 

2.9 

0.269 

— 

6.1 

0.3 

0.038 

0.024 

4.4 

0.8 

0.0 

2.78++ 

23.6++ 

12.3++ 

0.099+ 

688 

23.2 

14.5 

0.378 

2952 

16.1 

8.5 

0.125 

5120 

14.0 

9.6 

0.208 

3286 

10.7 

4.4++ 

0. 246++ 

2.25++ 

56.5++ 

42.0++ 

0.327 

2.45 

34.9 

30.9 

0.641 

— 

35.6 

24.6 

0.069 

0.309 

61.9 

43.3 

*  Amounts  In  litres  per  metre  of  shoreline  per  day. 

**  Amounts  in  grams  per  metre  of  shoreline  per  day. 

+  Amounts  in  milligrams  per  metre  of  shoreline  per  day. 

++  Estimated  from  established  flux  patterns  using  incomplete  data. 

♦  Possible  underestimation. 


Table  9 

Annual  Nitrogen  and  Phosphorus  Budgets  for  the  Lake 
Conway  System  for  the  Water  Year  1976* 


Nitrogen 

Phosphorus 

Inputs 

Airborne 

Combined  wet  and  dry  precipitation 

0.93 

0.125 

Surface 

Urban  stormwater 

0.84 

0.069 

Agricultural  runoff 

NS 

NS 

Undeveloped  land  runoff 

NS 

NS 

Subsurface 

Groundwater  inflow 

NS 

NS 

Seepage 

0.63 

0.024 

Septic  tanks** 

0.10 

0.006 

In  situ 

Nitrogen  fixation 

NS 

— 

Sediment  leaching 

(?) 

(?) 

Nutrient  recycling  (vascular  plants) 

(?) 

(?) 

Total  in 

2.50 

0.224 

Outputs 

Surface  outfalls 

0.22 

0.006 

Groundwater  recharge 

NS 

NS 

Fish  and  weed  removal 

NS 

NS 

Volatilization  and  evaporation 

(?) 

-- 

Denitrification 

(?) 

— 

Sed imentation 

2 . 28+ 

0.218+ 

Total  out 

2.50 

0.224 

Note:  NS  »  not  significant. 

*  Loadings  are  reported  in  grams  per 

square  metre  of 

lake  surface 

o  *  ~  —  —  o 

area  per  year  (g/nr-yr). 

**  Determined  by  calculation, 
t  Obtained  by  difference. 


i 


Table  10 

Monthly  Loadings  of  Nitrogen  and  Phosphorus  to  the  Lake 
Conway  System  for  the  Water  Year  1976* 


Date 

Storm¬ 

water 

Seepage 

Atmo¬ 
spheric  Outflow 

Storage 

Net 

Sedimentation* ** 

Nitrogen 

October  1975 

0.643 

0.487 

0.612 

0.00 

37.07 

-8.44 

November 

0.075 

0.067 

0.382 

0.00 

28.00 

+9.66 

December 

0.059 

0.052 

0.374 

0.00 

35.29 

-6.74 

January  1976 

0.042 

0.037 

0.366 

0.00 

11.86 

+23.94 

February 

0.094 

0.086 

0.392 

0.00 

18.65 

-6.156 

March 

0.296 

0.176 

0.442 

0.00 

23.01 

-3.48 

April 

0.347 

0.221 

0.467 

0.00 

22.14 

+1.87 

May 

1.386 

1.063 

0.928 

0.00 

29.49 

-4.11 

June 

1.138 

1.029 

0.904 

0.00 

30.77 

+1.85 

July 

0.906 

0.722 

0.742 

0.00 

26.38 

+6.72 

August 

0.480 

0.333 

0.538 

0.38 

42.68 

-14.61 

September 

0.773 

0.603 

0.676 

0.22 

26.05 

+18.80 

Phosphorus 

October  1975 

0.049 

0.018 

0.080 

0.00 

0.716 

-0.133 

November 

0.014 

0.002 

0.053 

0.00 

1.08 

-0.251 

December 

0.012 

0.002 

0.051 

0.00 

0.87 

+0.275 

January  1976 

0.009 

0.001 

0.051 

0.00 

1.09 

-0.159 

February 

0.013 

0.003 

0.054 

0.00 

1.13 

+0.030 

March 

0.021 

0.006 

0.060 

0.00 

1.04 

+0.177 

April 

0.026 

0.008 

0.063 

0.00 

1.06 

+0.077 

May 

0.101 

0.039 

0.118 

0.00 

0.60 

+0.718 

June 

0.092 

0.037 

0.115 

0.00 

0.75 

+0.094 

July 

0.071 

0.026 

0.096 

0.00 

0.38 

+0.287 

August 

0.030 

0.013 

0.070 

0.003 

0.37 

+0.122 

September 

0.060 

0.022 

0.088 

0.004 

0.47 

+0.070 

*  All  loadings  are  reported  In  grams  «  10  . 

**  Septic  tank  Inputs  were  treated  as  a  constant  inflow  for  these  calcula¬ 
tions.  Positive  sedimentation  represents  a  loss  from  storage  to  the 
sediments. 


APPENDIX  A:  STORMWATER  HYDROGRAPHS 

Stormwater  hydrographs  for  storms  were  measured  at  the  corner  of 
Gatlin  Avenue  and  Bumby  Drive,  Orlando,  Fla.,  on  the  following  dates: 

23  March  1978,  4  May  1978,  16  June  1978,  and  21  June  1978  (Figures  A1-A4, 
respectively) . 


APPENDIX  B:  COMPUTER  LISTINGS 


The  following  is  a  listing  of  the  hydrologic — materials  systems 
model  written  in  the  Continuous  Systems  Modeling  Program  (CSMP)  Language 
Also  included  is  the  Job  Control  Language  (JCL)  for  coupling  CSMP  to 
the  Gould  plotter,  using  a  program  in  the  Northeast  Regional  Data 
Center's  systems  library. 


B1 


( 


10*10*0)*'*****  BLANCH 

//arir tlB  00  O#»«A00ll3Ol «S1 -LOAOLla  orsp.*H» 

''  00  08N.GATOft,CSAPvio?.Dl8piaHR  3 

//PLUTLia  88  OaUIiAOOiiiSl'Sl^nlk.  iStHV*0i*P"*HK 

//sVsLifiBo8D0§s:s^8i  ?8^v!a^g?^2|firs^5MK 
^!!5?Sl2T  0?  oJShV™ 


BLANCHE*  ',CLAS8»M 


{^-8§yy|fpjgb;?gf?f5t.fe««fe.!.l 

y ^ T05F001  00  UNltSsf BDA# i^ACE*  (TRK ,  t> 

« _  PCB.irtECMiFB 


*0«), 


*|CL«*04} 

NE6L«i»04) 

NfcfL»a04) 


//8YSLINP  58  IN.0I!.P.C0LL1,P»») 

TiQ00  ONiT«SVijA/ £A  >11 

!%Issv“si?  jfh38s;sj!ile*f > 

;|Ppj;g8:sM-;?8ii?S«E?i5SIJ86:i?8J:}5f:J>St-a55i0) 

//srsiN  oo  * 


••♦•CONTINUOUS  SYSTEM  MQOtLlNG  PROGRAM**** 
***  VERSION  l.i  ••* 

TITLE  lake  CUNUAY  HYOROLOGIC-PHOSPHOKUS  model 
INI1IAL 


STATE  VARIABLE*  ***************** 


iSESS  .jS1!  3,488888! 

INCUN  bVlN  «  3650000}' 

•consianTs  used  In  my 


•COnsIaNTs  used  In  m?or 
CONST  rF«  1.0 

88884  u«»l*.8«‘*“  n 
8888!  331888  8  I18I888! 

CONST  MlVSUM  •  0.0 
CONST  AERqHT  ■  85,50 

mi\ 


In  MYftROLOGIC  EQUATION  tGMS  UNITS  EXCEPT  ZEREHT) 
1.5E  05 


CONST  A 
CONST  C 


1  * ” wvn  ■  v.w 

ICPOT  ■  0,11 
;PCO«  0.10 
JfiPMT  ■  0.92 


const  a2  ■ 

CONST  a3  ■ 


CONST  R3 
CONST  a  1 
CONST  ASE 
CONST  KR 

8888!  888 

CONST  AFL 


AR|Lal «0 

_ _  aP|HC  i  0.08 

CONST  kSfLOMbQ.1 
•  CONSTANTS  USED 
CONST  K3E0  ■  0,7 
A2  ■  0,2 
A 3  ■  0.2 


IN  MATERIALS  EQUATION  SMS  UNITS 


Jeep?  2lo.ose 

KRAlNl  ■  0,04 
RUNIN  »  0,39 

m :  km 


07 

130.0 


FUNCTION 


FUNCTION 


FUNCTIUN 


m 


FUNCtluN  BVAK  a  0, 

i.U: 


4*3 


0 1 4 • t  •  • 

ft* ,6 ■ Oft*** ( 


FUNCTIuN  AKOt 


iVA  6l?V 


mm* 

'S3,*A4  % 


81  » 


&!![»? 

Siflili: 

!:«:«: 

i:SS:ZS: 


\w 


i:? 

I:!* 


u: 

»! 

II; 

61. 


'  l|»# il*» l?i» •7*»^9a»  «!*  • 


FUNCTION  tKHT 


©iiiiaife 


;#s,}7;2ft;;sft.74,29' 

JM 


5**g£: 

. J: 

:i:iZ;SS: 
53:8:?!: 


f  ft 

ftftft 

•  M 

•  ft 
t  •  « 
ftftft 
ftftft 

ftftft 

ftftft 


•  • 


66,77* .  , 
8:8;:: 
8:8::: 


•***LAKE  VOLUME  EQUATION  ***»  * 


Q!N5SAfNi^Syft8^F»EEojlN*fiFLO*(  ♦  CPIN 
UOUT*EVAPtFLOUT*CPOUT*ORAIN 


(«**grounqnater  EQN  ***** 
6V*  2  JNTGRL  (*55lN  ^ 


IU :  mp-SMl • €IN  • 6uun 

G«lN"*RWWT*CVNl^*AR£A»5*KPkRC 

GFL0N--K6FL02* ( (KREL*V) *8V) 


«***  MATERIALS  EQN  ****  • 


QUTFLX  ■  ?CP0UT*VCl/VJ*tFL0Ur*yCi/V) 
SEQMNT  *  KSEO  *  VCI 

■ROtair . . 

NEN  N  ■  INF LUX* (PLANT I  ♦  8EDXN) 

jjS"  So rk|8H#iSPtWi;iW;- 

PCNTIN  ■CnENIN  /INFLUXJ*100. 


*KSEtPI 


SEOXN) 


I!m2nNU^LT*1.E*3»FIMT1M«63,,0UT0EL*1#O,PR0EL*1,O 

PSTSJff  SfS  JtI  , SEOXN. NENXN  ,  PCONC 
PREPaM  OELTMT,  height 

&TOP 


In  accordance  with  letter  froa  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsiaile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Blancher,  Eldon  C. 

Large-scale  operations  management  test  of  use  of  the 
white  a«ur  for  control  of  problem  aquatic  plants  : 

Report  2  :  First  year  poststocking  results  :  Volume  IV  : 
Nitrogen  and  phosphorus  dynamics  of  the  Lake  Conway 
ecosystem  :  Loading  budgets  and  a  dynamic  hydrologic 
phosphorus  model  /  by  Eldon  C.  Blancher,  II,  Charles  R. 
Fellows  (Department  of  Environmental  Engineering 
Sciences,  University  of  Florida).  —  Vicksburg,  Miss.  :  U.S 
Army  Engineer  Waterways  Experiment  Station  ;  Springfield, 
Va.  :  available  from  NTIS,  1982. 

59  p.  in  various  pagings  ;  ill.  ;  27  cm.  — 

(Technical  report  ;  A-78-2,  Rept.2,  Vol.4) 

Cover  title. 

"August  1982." 

"Prepared  for  U.S.  Army  Engineer  District,  Jacksonville 
and  Office,  Chief  of  Engineers,  U.S.  Army  under 
Contract  No.  DACW39-76-C-0076-002." 

"Monitored  by  Environmental  Laboratory,  U.S.  Army 
Engineer  Waterways  Experiment  Station." 

Bibliography:  p.  35. 


Blancher,  Eldon  C. 

Large-scale  operations  management  test  :  ...  1982. 

(Card  2) 

1.  Aquatic  biology.  2.  Aquatic  weeds.  3.  Conway, 

Lake  (Fla.)  4.  Weed  control--Biological  control. 

I.  Fellows,  Charles  R.  II.  United  States.  Army. 

Corps  of  Engineers.  Jacksonville  District.  III.  United 
States.  Army.  Corps  of  Engineers.  Office  of  the  Chief 
of  Engineers.  IV.  U.S.  Army  Engineer  Waterways 
Experiment  Station.  Environmental  Laboratory.  IV.  Title 
V.  Series:  Technical  report  (U.S.  Army  Engineer 
Waterways  Experiment  ition)  ;  A-78-2,  Rept.2, 

Vol-  *■ 

TA7.W34  no. A-78-2  Rept.2  Vol.4 
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